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Abstract 

We discuss a class of evolution equations equivalent to the simplest Universal Field 
Equation, the so-called Bateman equation, and show that all of them possess (at least) 
biHamiltonian structure. The first few conserved charges are calculated. 







1 Introduction 



Bateman's equation arose in a 1929 paper on hydrodynamics |jT| when he studied the prop- 
erties of the equation: 

^ = 'I^l4>tt - '^(PxMxt + <Pl4>xx = 0. (1.1) 

Here (p is a scalar field and subscripts denote partial derivatives. Many of the interesting 
properties of the Bateman equation are dealt with in ; its generalisations to higher dimen- 
sions are discussed in that and subsequent papers 0, ||, where they become the Universal 
Field Equations, a class of partial differential equations with intriguing Lagrangian formula- 
tions. Further analyis appears in |0. In the present case, we find that (|1 . 1|) has an infinite 
number of inequivalent Lagrangians. Any smooth, real-valued function homogeneous of 
degree one in the derivatives (p^ and (pt, and with arbitrary dependence on (p, works as a 
Lagrangian for A. 

This is straightforward to prove. Suppose C{(p, (p^, (pt) is a candidate Lagrangian. By 
Euler's theorem on homogeneous functions: 



1^ + = (1.2) 



It follows easily that: 



"dcpl ^dcp^dcpt 
'd(p^d(pt ^ dcpl 



Now, if we expand the Euler-Lagrange expression. 



and use the homogeneity properties ( pT3| ) the equation of motion becomes: 

1 d^C 

= ~~r a- {4>l4>tt - '2(px(pt(pxt + (Pt(pxx) ■ (1-5) 

(px(ptO<pxO(pt ^ ' 



This is equivalent to ( |1 . 1| ) on extremals of L. 

Note that the homogeneous equations in (|1.3| ) imply that L is always singular. Also, 
notice that if there is no dependence in the Lagrangian, the Euler-Lagrange equation is 
in the form of a conserved current, and so there are an infinite number of conservation laws 
for ( |1 ■ 1|) ■ We will use this observation to construct Hamiltonian conservation laws that are 
properly in involution. 

It should be noted that much is already known about the Bateman equation from the 
classical theory of partial differential equations. For example, in common with other quasi- 
linear second-order equations 0, A can be linearised by a Legendre transform [Q. With the 
change of variables, 

^ = 0i 

u{i,r]) = xi\yr\- (1.6) 



Lll) is equivalent to: 



+ 2^770;^^ + r/^cj^^. (1.7) 

This linearised equation supplies an implicit solution for A = 0. According to a theorem of 
Kumei and Bluman [Q, the linearisability property is due to the existence of a first-order 
generalised symmetry of (|1 . 1| ) whose characteristic Q depends on an arbitrary solution of 
the linear equation: 



From this we deduce that any Q homogeneous of degree zero or one in the first derivatives 
of is a characteristic for a generalised symmetry of ( |1.1D . 

The task at hand is to study a more formal matter, namely whether there is a Hamiltonian 
structure or structures associated with the various Lagrangian formulations of the Bateman 
equation, and how its integrability features in such a formalism. The basic reference for the 
techniques used is M. 
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2 Hamiltonian Analysis 



The Hamiltonian and multi-Hamiltonian descriptions of the Born-Infeld equation have al- 



ready been studied by Nutku and various collaborators |jTO|, |TT| and the aim here is to perform 
a similar analysis of the Bateman equation. The results of these papers can be adapted to 
the present case simply by using the A = limit of the (unique) Born-Infeld Lagrangian, 



CBi = ^\ + ct>l-^l (2.1) 



This essentially repeats the work in |TT| but, as we have seen, the set of admissible La- 
grangians for (|1 . 1|) is far larger than this. We can discuss Hamiltonian formulations corre- 
sponding to a large class of these non-standard Bateman Lagrangians and in general we will 
find behaviour quite different from the Born-Infeld case. 

We need to find an evolution equation (or system of evolution equations) of the form, 

Ut = J5rL[Ml (2.2) 

where ?i[u] is the Hamiltonian functional of the dependent variables u,v and is a skew- 
adjoint differential operator. In this context, "skew-adjoint" means that 

/ A.JBdx = [ B.J*Adx = - f B.JAdx, (2.3) 
Jn Jn Jn 

assuming that the support of A and B on the region Q allows formal integration by parts. 

The associated Poisson bracket is: 

{V, Q} = J 6V.J6Qdx. (2.4) 
This bracket must satisfy the Jacobi identity, which can be verified using Olver's convenient 



method using functional tri- vectors ||12|, We require that this system is equivalent to (|1 . 1| ) 
under some suitable change of variables. 

Ideally, we would like to be able to find two such structures for our evolution equations. 



Then, provided some compatibility criteria are satisfied, Magri's theorem |]T3 guarantees the 



complete integrability of the equation and allows construction of the Lax pair representation 
and the associated conservation laws. 
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The most direct approach is suggested by the fact mentioned earher, that if the Bateman 
Lagrangian has no (j) dependence, ( |1.1|) can be expressed as a conservation law: 

The biHamiltonian stucture of equations expressible in the form of a conservation law was 
demonstrated by Nutku As mentioned above, there are an infinite number of ways of 



expressing the Bateman equation as a divergence-free current, and there is a biHamiltonian 
structure associated with each. 

Following [0], the divergence (|2.5| ) may be rewritten as the smoothness condition for a 
function i/j. 

A = (2.6) 

In accordance with the standard Hamiltonian technique, the field (f) is mapped to a pair of 
variables u, where 

u = <px, 

V = i)^. (2.7) 
The t-derivatives can then be expressed in terms of these new variables: 

0i = U{u,v), 

A = Viu,v). (2.8) 

Consequently, the equation of motion and the smoothness condition for (p can be summarised 
by the condition of closure for the pair of exact 1-forms, 

d(f) = u dx + U dt, 

dip = v dx + V dt, (2.9) 
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namely that: 

ut - D., (U) = 

Vt-D,{V) = 0. (2.10) 

Given our assumptions about the nature of the Lagrangian, we can deduce the form of 
U and V. From ( p. .31) , we know that the derivatives of £ with respect to (f)^ and (j)t are 
homogeneous of degree zero. So, let us write ( p.6|) as, 



(2.U) 



where F and G are smooth, invertible functions. Solving the first of these relations for (pt 
tells us that, 

U{u,v) =u F-\v), (2.12) 

and substitution into the second, 

V{u,v) = G(^F-\v)). (2.13) 

At this point, we attempt to impose the Hamiltonian structure. A nice result about this 
type of equation |T3] is that it admits the canonical, "fiat" Poisson structure defined by the 



structure matrix, 

J, = ( t ) ■ P.14) 

and a Hamiltonian Hi dependent only on u and v, subject to the conditions, 

which in turn impose compatibility criteria on U and V. 



In this example, the compatibility requires that: 

d_ 

dv 



G'{f~\v))^{f-\v))=F~\v). (2.16) 
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Bearing this in mind, we end up with a general Hamihonian framework for the Bateman 
equation. The Hamilton equations are easily seen to reduce to the form, 

ut = {uF'\v)) , 

vt = F-\v) V,. (2.17) 
The structure matrix is J7i and the Hamiltonian is: 

Hi = u G (^F-\v)) . (2.18) 
There is also a biHamiltonian structure associated with these equations. For this "con- 



served current" type of system there is always a conserved quantity of the form |14], 



Ho = uv. (2.19) 

Following the standard program for biHamiltonian systems we postulate this as a con- 
served Hamiltonian which reproduces the equations of motion when used with a second 
structure matrix J'2: 

JiSHi = J2Sno. (2.20) 

Provided this second structure is "compatible" with the first, meaning that any linear com- 
bination of the two is also a Poisson structure, this construction will give rise to a (hopefully 
infinite) sequence of conserved quantities Hk such that, 

JiSHk = JiSUk-u (2.21) 

which can be generated by the recursion operator: 

n = J^J^y (2.22) 

We know that the general form of J^-, as the most general first-order, skew-adjoint, 
matrix differential operator is, 

^ / ImDj, + 2pDa; + (p + g)^ \ ,^ 

^2 1^ 2pD, + (p - q\ 2nD, + j ' ^^■^''> 
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where m, n,p, q are functions of u and v. By using the condition ( p.2(]| ) and demanding that 
it reproduces ( p.l7|) , restrictions on the form of J'2 can be found. Nutku [|l4l proves that 
these restrictions also ensure the closure of the Jacobi identity. 

It is difficult to write down a convenient general form for the most general structure 
so derived, but a useful simplifying assumption is to take n = 0. Then we find a possible 
candidate for the second structure given by, 



k 



m 



n = 0, 



P 



q = \F-\v) + k\ 



(2.24) 



where k and k' are constants. It is straightforward to ascertain that this structure is compat- 
ible with the Ji. the combination Ji + J2 is skew adjoint and satisfies the Jacobi identity. 
So J2 can be used to find a second conserved Hamiltonian and thence a whole hierarchy of 
conserved quantities. 

The associated recursion operator is the pseudodifferential operator. 



7^ = J2J{ 



-1 



\^ F-\v) 

using the values of p, q, and m in ( p.24|) . 

It is known H that a system of evolution equations of the form. 



(2.25) 



ut = K[u], 



(2.26) 



which admits a recursion operator, admits a Lax-type representation. 



U = [A, L] 



(2.27) 



if the L is the recursion operator ( p.25|) and A is the Frechet derivative: 



^=Dk= -K[uJ 



£=0 



(2.28) 



Here, denotes a one-parameter (e) family of perturbations of u. 
In our example, the Frechet derivative is: 



A=( (^-H^)-) H^-\^)-j 1 ^2 291 



The argument of the A should be inserted as indicated by the dots. Substitution of ( p.29|) 
and ( |2.25| ) into ( p.27| ) reproduces the Hamilton equations ( |2.17| ). 



Now, using the recurrence ( p.21| ), we can start to calculate conserved quantities. The 
first conserved Hamiltonian is found by considering: 



It turns out to be: 



J2Sni = Ji6n2. (2.30) 



H2 = f {F-\w)f + '^^G{F-\w))^ dw. (2.31) 

Continuing the process, we find the third charge: 

^3 = f \f-\w) [u{F'\w)f + '^^G{F-\w)) + ^ f{F~\w')fdw'^ | dw. (2.32) 

The subsequent members of the hierarchy become rapidly more complicated. 

It is worth noting that the conserved-current type of equation found by Nutku admits 
a third independent Hamiltonian structure, compatible with the first two, with 7^2 as the 
relevant Hamiltonian density. Also, while the sequence of conserved charges in principle 
continues ad infinitum, there may be certain particular cases when the sequence terminates 



or repeats itself after a finite number of steps. See [11] for an example of this behaviour 
arising in the biHamiltonian hierarchies of the Born-Infeld equation. Essentially the same 
results can be found for the Bateman case by choosing the Lagrangian ( p.l|) with A = 0. 
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